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Two-loop computation of a finite volume running coupling on the lattice 
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'^Institut fiir Physik, Humboldt Universitat 
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In pure SU(2) gauge theory we compute the two- loop coefficient in the relation between the lattice bare coupling 
and the running coupling defined through the Schrodinger functional. This result is required to relate the latter 
to the MS-coupling in our programme to compute as. In addition it allows us to implement 0(a) improvement of 
the Schrodinger functional to two-loop order. The two-loop /3-function is verified in a perturbative computation 
on the lattice, and the behavior of expansions in the standard and in the Parisi-improved bare couplings are 
investigated beyond one loop. 



1. CONTEXT OF THE CALCULATION 

In this contribution we report on a perturba- 
tive 2-loop computation of the running coupling 
constant 5^(i), defined through the Schrodinger 
functional M , in terms of the lattice bare coupling 
go- Such a result is required as part of our project 
of computing as for QCD. The strategy there is 
to first numerically construct the continuum limit 
oig for some large values of L~^ in physical units. 
These are then connected vifith the more conven- 
tional 5j^(q) by renormalized perturbation the- 
ory. This is expected to be a well-behaved appli- 
cation of the series, as the scales L~^ and q are of 
the same order and can be chosen well above the 
nonperturbative scales. With growing statistical 
accuracy, especially on parallel computers, the 2- 
loop coefficient in this conversion is required to 
achieve a reasonable balance and control of the 
systematic errors from truncating the series. The 
connection between g and gj^ is established by 
expressing both of them in terms of the lattice 
bare coupling go. In this way, apart from avoid- 
ing problems of a direct connection via dimen- 
sional regularization in the presence of a finite 
background field, we have the additional benefit 
of getting the two- loop 0(a) improvement coeffi- 
cient to speed up the continuum limit ||^. Here 
we present the results of a first such computation 
of 5 for SU(2) pure gauge theory Q. For the con- 
nection with gj^ see Peter Weisz's contribution 
to these proceedings. 



2. FINITE VOLUME RUNNING COU- 
PLING 

The definition of 'g from the Schrodinger func- 
tional is discussed in more detail in ||l|,|[, and we 
are brief here. The effective action F is given by 
the path integral over SU(2) gauge fields U{x,ii) 
on a lattice with extension L in all four directions, 
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While space is periodic, the fixed boundary values 
in the time direction, U{x, A:)|j.o^q = explrjT^/iL} 
and U{x,k)\xO^]^ = exp{(7r — ri)T^/iL\, induce 
an ?7-dependent abelian background field. The 
action S is the standard Wilson action modified 
for the surfaces at a;° — 0, L, 
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^[f^] = -E^(pMi-^(p)}' 
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where the sum runs over all oriented plaquettes p. 
The weight w{p) is unity except for the the electri- 
cal plaquettes touching the fixed field boundary, 
where it equals 



c,(5o) = l + c(%^ + cf'.go' 



(3) 



Its tuning allows for the suppression of 0(a) arti- 
facts arising from the presence of the boundaries. 
Following Symanzik, we perturbatively determine 
the coefficients in 
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ct{go) = l + cY>gt,+c^i>gl 
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The response 
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defines g, which runs with the finite system size L 
as the only scale in the problem. Here ' denotes 77- 
derivatives, and Fq is the exactly known classical 
limit (50 -^ 0) of r. 

3. PERTURBATIVE CALCULATION 

A suitable gauge fixing with the introduction of 
ghosts and a systematic expansion of the fiuctu- 
ations of the [/-field around the induced back- 
ground lead to a well-defined lattice perturba- 
tion expansion [|l|. In particular, all modes are 
quadratically damped with the present boundary 
conditions. In the resulting expansion 

f{L) = 5o' + mi{L/a)g^o + ^^iLjo^gl + • • • , (6) 
the coefficient m\ was given in R] and, at 2 loop 
order, TO2 is reported here and in more detail in 
1^ . Evaluation of (the ry-dependence of) the order 
of 10 vacuum 2-loop diagrams is required for vfii. 
The computation is involved, because all prop- 
agators have to be computed numerically with 
the finite background field in place, and because 
the 3- and 4-gluon vertices have a large num- 
ber (~ 100) of terms. Hence the diagrams are 
summed numerically for Lja < 32. To get to this 
size on workstations in a reasonable time, all sym- 
metries were used eventually to reduce the num- 
ber of terms. On smaller lattices they were how- 
ever checked to hold first. As further error checks 
independence of a continuous gauge fixing param- 
eter was verified, and the two authors of Q car- 
ried out independent calculations for the smaller 
lattices before optimizing one of the codes. 

The resulting column of values of m2 vs. L/a 
was analyzed as described in B . The asymptotic 
result is 
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2ba\n{L/ a) + 0.20235 + 0(0^/1^) (7) 



26iln(L/a) + 0.01607 +0(aVi^)- (8) 

Here bo, 61 are the universal coefficients of the /?- 

function 
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bo ^777^, bi^T^^. (9) 



The value of 61 was confirmed in our analysis of 
m2{L/a) to about 1 part in 10^, and to our knowl- 
edge this is the first such check on the lattice. 
The lattice artifacts of 0(a) have been canceled 
in mi , TO2 . This requirement fixes the improve- 
ment coefficients 






-0.0543(5), a 



(2) 



-0.0115(5). 
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4. EXPANSION IN THE BARE COU- 
PLING AT TWO LOOPS 

As mentioned, the result of the previous sec- 
tion is needed to connect g and gj^ via gg. As 
a byproduct, we have the opportunity to study 
the quality of the lattice perturbation expansion 
in go for a physical quantity. Eqs. (||), (|^), (^) 
can be used to relate g at a scale a/s, s = 0(1), 
in the continuum with the bare lattice coupling 
go associated with the cutoff a in the regulariza- 
tion in use. On the one hand, no large logarithms 
are expected here as scales of the same order are 
connected. On the other hand, the expansion in- 
volves non-universal lattice quantities and can- 
not be systematically improved for fixed physical 
scale by approaching the continuum limit. The 
precise derivation of the relation, although sim- 
ple in result, is shghtly tricky to argue Q: One 
first uses (|6|) in a regime close to the continuum 
limit, where go, goln(L/a) are small in the ex- 
pansion, and d^ jl? is negligible. The resulting 
g^(L) is re-expanded in terms oi'g^{a/ s) by using 
continuum perturbation theory. This procedure 
is formally equivalent to dropping powers of a/ L 
in (ra), (@) and then putting a/L = s under the 
logarithms. 

In terms of a{q) = g'^ (L) / (Att) , L = q' 
ao — go/{4:TT) one gets for s = 1 

a{a^^) =ao + 2.543 ag 



and 
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The rather large coefficients here are consistent 
with the experience, that go used in this way is 
not a good expansion parameter. An alternative 
way of using the series is to fix s such that there 
is no 1-loop term, and this leads to 

a(8.83 a"i) = ao + 1.287(l)ai] + 0{a^). (12) 

Note that this choice simultaneously produces a 
small 2-loop coefficient, which is nontrivial and 



leaves the chance, that the series may be better 
behaved with this rather large scale shift s — 8.83. 
Recently the use of modified bare couplings has 
become popular to try to amend the bare series 
J6| . A simple proposal is 



Qo = ao/P, 



(13) 



where P is the average plaquette in infinite vol- 
ume. Using the known series for P, one finds 
analogously to (|2|) 



Q(1.17a-i) = ao + 0.951(1)5^ + O(a^) 



(14) 



Here the 2-loop coefficient is reasonable also, and 
the scale factor is much closer to one. 

We now want to compare the expansions (12) 
and dl4) with nonperturbative results for "g. We 
consider /3 = 2.85, where ao and (5o are available, 
and the scale a^gg is about 8 GeV if the SU(2) 
potential is matched with nature H . The expan- 
sions produce g at 1.17 a^g5 and 8.83 a^g5, which, 
for the comparison, we both evolve to 10a2^g5 
with the numerically controlled /3-function of g 
and negligible error. At this scale g^ itself is 
known numerically , and all values are collected 
in table 0. The 2-loop results are only off by about 
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Estimates for 


a(10 


^2.85) 


a 




method 


0.1135(8) 




nonperturbative 


0.1098 




1-loop in QfQ 


0.1115 




2-loop in ofQ 


0.1110 




1-loop in ao 


0.1128 




2-loop in ao 



experience depends on the renormalized quantity 
(here g) that is computed and need not be uni- 
versally true. 
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1 and 2.5 error margins for ao and ao respectively. 
When using the 1-loop term to fix the scale for 
which to apply bare perturbation theory, we find 
both a reasonably sized 2-loop term and approxi- 
mate agreement with the "exact" result. Perhaps 
a bit surprisingly, this works for the ordinary bare 
lattice coupling in a qualitatively similar fashion 
as for ao, with just the scale factor being rather 
large. One has to be cautious, however, that this 



